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SUMMARY 


A relatively short method of calculating the lateral motions of an 
airplane due to random atmospheric turbulence is presented. The gust 
velocities are represented as equivalent rigid-body rotations of the 
airplanes; namely, rolling gusts, yawing gusts, and side gusts. Random 
distributions of gust velocities across the span are taken into account 
in defining the rolling and yawing gusts. Complex stability derivatives 
are used to account for the random distribution of side gusts along the 
fuselage and vertical tail and the lag effect incurred as the airplane 
penetrates the gusts. The suggested gust spectrum is based on a simple 
analytical expression which fits available measurements of atmospheric 
turbulence. A 45-step sample calculation procedure for obtaining the 
response of the airplane in each degree of freedom is presented in 
tabular form. 


INTRODUCTION 


Most calculations of the response of airplanes to gusts have been 
made on the assumption that the effect of the gust on the airplane is 
approximately equivalent to the effect of a rigid-body motion of the 
airplane producing a distribution of angle of attack similar to that 
caused by the gust. On the basis of this assumption, disturbances in 
the form of rolling gusts, yawing gusts, and side gusts have been employed 
in calculating lateral response to gust disturbances. This approach is 
convenient because the standard aerodynamic stability derivatives which 
are used in airplane stability calculations may also be used to determine 
moments caused by the gust velocities. This approach as usually applied 
neglects effects due to lag in penetration of the gusts by different parts 
of the airplane and, because linear gradients of the gust velocity along 
the span are assumed, it cannot account for the random spanwise distribu- 
tion of gust velocities encountered in flight through atmospheric turbu- 
lence. Furthermore, the relations between the magnitudes of the rolling, 
yawing, and side gusts required to produce effects similar to actual 
atmospheric turbulence are not known beforehand. 
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In reference 1, a theoretical method for calculating the lateral 
response of an airplane to atmospheric turbulence has been proposed 
which account's in a rather complete manner for the effects neglected or 
approximated in previous methods . This method uses an approach somewhat 
different from that described in the preceding paragraph in that the 
forces and moments applied to the airplane by gusts are determined in 
power -spectral form in terms of the horizontal, vertical, and side com- 
ponents of gust velocity. 

In the present report, it is shown that the more conventional 
method of assuming gust-velocity distributions equivalent to rigid-body 
motions of the airplane may be refined to provide results equivalent to 
those given by the method of reference 1. The refinements consist in 
replacing some of the constant aerodynamic stability derivatives with 
complex quantities to account for the gust penetration effects and 
determining the correct relations between the spectra of rolling, yawing, 
and side gusts to yield results in agreement with the more exact analysis . 
The present method requires somewhat simpler calculations than the method 
of reference 1 and provides a clearer physical picture of the relations 
between the various sources of lateral gust disturbances . 


SYMBOIS 


b 


wing span 


D 


nondimens ional operator, — 


U dt 


G 


matrix containing stability derivatives relating air- 
plane moments and forces to gust velocities 


G 


alternate form of G containing frequency-dependent 
stability derivatives 


h 


height of center of pressure of vertical tail above 
X-axis of airplane 


I 


imaginary part of a complex number 




nondlmensional radius of gyration about X-axis 


nondime ns ional radius of gyration about Z-axis 



nondimensional product of inertia 
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k' = aii/u 


L 

integral scale of turbulence 

l t 

tail length between airplane center of gravity and 

quarter -chord point of mean aerodynamic chord of ver- 
tical tail 

P 

rolling velocity, d0/dt 

<1 

1 2 

dynamic pressure, ^pU 

r 

yawing velocity, d^/dt 

s 

wing area 

T 

time interval over which power spectrum is evaluated 
(eq. 18 ) 

t 

time 

U 

relative velocity between airplane and general air mass 

u 

velocity along X-axis 

V 

velocity along Y-axis 

w 

velocity along Z-axis 

X,Y,Z 

three orthogonal reference axes of airplane 

c L 

lift coefficient, ^ift 

qS 

c 1 

. . Rolling moment 

rolling -moment coefficient, — 

qSb 

c n 

. . . . Yawing moment 

yawing -moment coefficient, — 

C Y 

side -force coefficient. Side force 

qS 

a 

angle of attack 

3 

angle of sideslip 


k 
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3 ' = b/L 

7 flight -path angle 

A matrix containing airplane equations of motion in still 

air 


A 

P 

0 

0 

O) 

a ) 1 = a)b/U 


wavelength, 

O) 

airplane relative density factor, 


Mass 

pSb 


density of atmosphere 
sidewash angle 
power spectral density 
angle of roll 
angle of yaw 
circular frequency 


Stability derivatives of airplane are indicated by subscript nota- 
tions j for example , 



' y p = aT 


Subscripts : 

o general air mass 

g gusts 

W wing 

F fuselage 


T 


vertical tail 
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Matrix notation: 

determinant or absolute value of quantity 
rectangular matrix 
column matrix 
row matrix 


[] 
{ } 


L J 


Bar over a quantity denotes mean value. Asterisk denotes complex 
conjugate . 


THEORY 


Equations of Motion 

The equations of lateral motion for an airplane (ref. 2) are 
given by 


2hK x 2 D 2 0 o - |Cj p D0 - 2^K XZ D 2 * 0 - ^Djr - C^p = 0 

-2vKwP% - + 2hK z 2 D^ 0 - - C^p = 0 

-iCyp^0 “ ^L^o + 2(j.E>4r 0 - TjCy^Ity ” Cl tan y \|r 0 + 2 ^DP q - Cy^p = 0 


( 1 ) 


where the subscript o appearing in the inertial and weight terms is 
used to denote angular displacement with respect to an absolute system 
of axes fixed in the general air mass. In calculations of the motion 
of an airplane in still air, the angular displacements and velocities 
appearing in the aerodynamic terms are identical with these values. 

When flying in turbulent air, however, the airplane is subjected to the 
motion of local air masses, generally referred to as gusts. The rela- 
tive linear and angular velocities of the airplane with respect to the 
local air mass each may be considered as made up of two parts: 
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D0 = D0 O + D0 g 

DJr = D|r 0 + Dlr g ► 

P = P 0 + 3 g 


( 2 ) 


where "the subscript g is used, to denote gust velocities . Substituting 
equations (2) into equations (l) and transposing the terms resulting from 
gust disturbances to the right-hand side of the equation gives the result 
written in the convenient matrix form: 



to 


'^g 

W- 

*0 

-[«]■ 

Ofg- 


0o 

V i 


0g 

\ . 


The matrix 


( 5 ) 



2^K x 2 D 2 - |c Zp D 

-2^xZ d2 ' 

- ic 7 D 
2 6 r 


H- 

- 2 M iC XZ I)2 “ p-'np 0 

2hK z 2 d2 - 

ic n D 
2 n r 

-s 


-k: Y D - C T 
2 Y P L 

1 

P 

1 

tan 7 2)aD - Cy 


is the familiar "still air" rigid-airframe characteristic equation, 
the matrix 


00 


and 


al- 


¥ l V 


iC 

2% 


ic Y 

2 Y p 


¥ l r 

% 

2° n r 

Cn P 

2° Y r 

% 


( 5 ) 
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gives the relationship between the aerodynamic moments and forces 
resulting from gust velocities encountered. 

In the classical method of treating simple Isolated gust inputs 
(for example, ref. 3 ), the elements of the G matrix would be the 
"still air" stability derivatives, the values in the gust-velocity 
matrix ( D 0 g , D|fg, and would be those appropriate to step or ramp 

functions, and the solution would be defined by 


0 o (o>) 


D 0 g (a>) 

< t Q (co) 

> = [am ]" 1 [g] < 

Dt|fg(o)) > 

!0 o M 




where the inverse of [a"| with D = ij^o). Unique time 

responses could then be obtained by taking the inverse Fourier trans- 
form of each solution. 

However, when the responses of the airplane to continuous random 
gusts are to be considered, the gust velocities can be defined only in 
a statistical (power-spectral) manner and the resulting moments and 
forces will in turn be related only in a statistical sense. This means 
that the elements of matrix G cannot be evaluated in the usual sense 
and the effect of random distributions of gust velocities along the 
fuselage and across the wing span must be taken into account. 


Forces and Moments Due to Turbulence 

In the application of the method to random turbulence, the sources 
of the forces and the moments on the airplane must be considered. As 
shown in reference 4, yawing and rolling moments on the wing result from 
gradients of the horizontal and vertical gusts . Wing rolling moments 
and moments and forces on the fuselage and vertical tail are produced 
by side gusts. Although the wing moments cannot be determined at any 
instant as a function of the gust velocities measured at the center 
of gravity, the power spectra of these moments have been determined 
in reference 4 as a function of the power spectrum of the vertical gust 
velocity as measured at one point. In isotropic turbulence, the 
spectrum of the side gusts measured at a point is identical to that of 
the vertical gusts. For this reason, the spectrum of the yawing and 
rolling moments may be related equally well to that of the side gusts. 
This procedure is used in the present analysis. 
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Spanwise gradient of vertical gust .- In general, the rolling moment 
acting at any instant results from a random distribution of vertical- 
gust velocity across the span which has some average linear spanwise 
gradient. An effective gradient which produces the rolling moment due 
to this random spanwise distribution may be defined by the relation 




“ c z( w g) 

w 


which is derived from the equation 


c *<"e> * |(°! P ) W ^ 


(7) 


where Ci(v s ) is the rolling moment due to vertical gusts, (^i ) i s 
the wing damping -in-roll stability derivative, and D 0 g is the equiva- 


lent rolling -gust gradient. Likewise, the vertical-gust distribution at 
any given instant also produces a yawing moment. In reference 4, this 
moment was assumed to be in phase with the rolling moment and was deter- 
mined from the formula 


Cn( w g) 



( 8 ) 


Substituting equation ( 7 ) into equation ( 8 ) gives the yawing moment in 
terms of the effective gust gradient: 



(9) 


Spanwise gradient of horizontal gust .- Similar arguments may be 
used to show that a random distribution of horizontal gusts across the 
span of the wing produces both rolling and yawing moments which are 
assumed to be in phase and related by 

/ c n r \ 

c n( u g) - ( c? ) C l( u g) 

V 1 r /w 


(10) 
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In terms of the equivalent yawing -gust gradient Ityg, these moments are 
given by 


c K u g) " a( c *r) w ^ 

(ID 

c n( u g) ■ |( Cn r)u n| '8 

(12) 


These formulas (eqs. ( 7 ) to (12)) are in a convenient form for 
application to the present analysis in which the gusts are considered 
as equivalent rigid-body motions of the airplane . Moreover , since only 
the wing moments axe influenced by the horizontal- and vertical -gust 
distributions across the span of the wing, these formulas completely 
account for the effect of these gusts on the lateral motion. 


Side gusts . - The only remaining source of gust disturbance is the 
side gust. In the present report this disturbance is expressed as an 
equivalent sideslip 0 g = v g /U, but the method of calculating its effect 


is essentially the same as that used in reference 1. Effects due to 
the difference in time at which a given gust encounters different sec- 
tions of the airplane, known as gust-penetration effects, are accounted 
for in reference 1 by considering the relations between the aerodynamic 
force and moment coefficients and the gust inputs measured at the center 
of gravity as frequency-dependent transfer functions . These relations 
are herein converted to frequency-dependent stability derivatives Cy^ , 


and C 1 


the flight velocity U. 
are given in appendix A 


by multiplying the expressions used in reference 1 by 
Expressions for these stability derivatives 


To an observer in the airplane the side -gust disturbance appears 
as a random velocity distribution which moves along the fuselage and 
vertical tail at the mean velocity U. If this speed in terms of body 
lengths per unit time is large (as is usually the case with most air- 
planes), these gust velocities will not change appreciably during their 
"time of exposure" to the fuselage. When the condition is satisfied, 
the forces and moments acting on the fuselage are uniquely, rather than 
statistically, related to the side -gust distribution along the flight 
path, and phase relations can be correctly accounted for by treating 
the frequency -dependent derivatives as complex quantities. 
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Matrix Solution to Equations of Motion 

The effects of the gusts on the wing and the fuselage -tail combi- 
nation may be incorporated into a matrix similar to the simplified 

G matrix of equation (5). If this matrix is denoted by G(cd) and if 
only those force and moment coefficients are included which will have 
a significant effect (for example, side forces due to rolling and 
yawing are usually negligible compared with side force due to sideslip), 
the modified G matrix is defined as 





[V" , j 

Lt 

[S(«.)j - 

!(%)„ 

!( Cn r) w 

[%(“>>] 

FT 


0 

0 


l FT 


In terms of this matrix, the solution to the equations of motion given 
by equation ( 3 ) is obtained formally by inversion: 




D0 g (a>) 

V 0 (o>) 

' = £a(cd)J _1 [S(aj)j < 

D*t r g(a^)> 

* ✓ 


y - 1 , 


(14) 


The inverse of the A matrix given by equation (4) is shown in 
appendix B to consist of transfer functions relating the response in 
lateral angular displacements to a sinusoidal rolling moment, yawing 
moment, or side force of unit amplitude: 



£-(<») 

/(a)) 

/(a)) 

c l 

C n 

C Y 

Uo) 


/(a)) 

Cl 

c n 

C Y 


Wi 

—(a)) 

C Y 


( 15 ) 
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Inasmuch as the random nature of the moments resulting from random 
turbulence mait.es it necessary to place the result given by equation (l4) 
in a power-spectral form, both sides of the equation must be squared. 

If equation (l4) is squared in its present form, cross-power terms 
between all the forces and moments will appear in the final result. 

However, if the product of [a ] -1 and [~G(cd)J is taken beforehand, 
then when the product matrix is squared these cross -power relationships 
do not appear explicitly since their equivalent effect has been taken 
into account in the multiplication process . It is shown in appendix C 
that the result obtained by first multiplying these two matrices before 
squaring is the same as the result obtained by squaring each matrix 

separately and including cross -power terms. 

Each element in the matrix product of [aJ - ^ and [~G(co)J has the 

form of a transfer function relating one of the airplane response 
quantities to one of the gust components. This relationship is indi- 
cated as follows : 




i- (ox) 
^g 

1(0,) 

Pg 

[a-Igm] = 



f(»> 



(u>) 

P g 


(16) 


In terms of this matrix product, the relations between the power spectra 
of the airplane response and the gust components are given in the form 


H 


1 

1 

A. 

' 5 ^ 


> = 


* 

% 

V 



3 


where the power spectrum of a 


0 

^g 

2 

1 

0 

p g 

2 


X 


2 

Pg 

2 

< 

X* 

3 

2 

3 

P g 

2 


$ 3g 


quantity X is defined by 


( 17 ) 
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lim 
T -» 00 



( 18 ) 


The assumption has been made that any cross power between the three gust 
components is either zero or negligible. Justification for this assump- 
tion is made in a subsequent section. 


Finally, the gust inputs in isotropic turbulence are specified in 
terms of a single quantity. The spectrum of 3g is selected for this 
purpose inasmuch as Pg is directly related to the linear gust compo- 
nent Vg for which the spectrum is available from turbulence theory. 
Therefore, the final form of the equations is given by 


V 


~ 0 

2 

0 

2 

0 

2 



2 






p g 



p g 


v 

> = 

* 

2 

♦ 

^g 

2 

t 

0g 

2 

< 

^g 

0g 

2 

V 

> 


1 

2 

0 

2 

p 

p g 

2 


1 

i 


( 19 ) 


Gust Spectra and Their Relationships 

The following relationships for the gust spectra and gust velocities 
are based on the assumptions of homogeneous isotropic turbulence : 


(<d) = ® v M = U 2 ® p (a>) 
8 8 8 



(20) 


Physically, these relations state that, regardless of the motions or the 
direction of travel of the airplane, the vertical- and side-gust com- 
ponents measured at the same point on the airplane have the same spectrum 
and that the mean-square value of all three gust components is the same . 

For the purposes of equation (19), the power spectra of gust gradi- 
ents as a function of the gust spectra measured at one point are required 
in the form of the ratios of the power spectra of D 0 and EH|r to the 
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power spectrum of Pg when the latter is known. Taking the square of 

the absolute value of equation (7) and. considering the equation to apply 
at any given frequency yields 





( 21 ) 


Dividing both sides of equations (21) by the power spectrum of the side 
gust or vertical gust (see eqs . (20)) as measured at the center of grav- 
ity gives 


®D0 


g 


4U 2 * c l(»g) 


(%); 


*w 


g 


At any given frequency, then 


D0 


g 


4U 2 %(w g ) 


(%)w 


( 22 ) 


w 


g 


Expressions for both power spectra on the right-hand side of equa- 
tion (22) are available in the literature. An expression based on 
measurements of turbulence in wind tunnels (ref. 5) which appears to 
fit well most of the available data from measurements of atmospheric 
turbulence (for example, refs. 6 and 7) is given in terms of Pg by 



g 



v g 2li 1 + 3k’ 2 
rtU (1 + k ' 2 ) 2 


(23) 


where k 1 = cjuL/U and L is the so-called scale of turbulence. A plot 
of equation ( 23 ) is given in figure 1 . 

Calculations based on analytical expressions (similar to eq. ( 23 )) 
for the gust spectra measured at one point are given in reference k for 
the power spectra of the coefficients of the rolling and yawing moments 
on wings of arbitrary span which are subject to continuous isotropic 
turbulence. These spectra, which take into account the random distri- 
butions of gusts across the span and along the flight path, are given 
in reference k for various values of 0 * = b/L and for four spanwise 
lift distributions on the wing. Since the effect of different lift dis- 
tributions was small, only one distribution (the rectangular distribution) 
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is considered and the power spectrum of rolling -moment coefficient due 
to vertical gusts on the wing is denoted here by 


Wg) 


itu3 


x Quantity plotted in fig. 7( a )> ref. 4 


The quantity plotted in figure 7( a ) of reference 4 is 




W B 2LC !p7 U5 ' 

Dividing the above equation by equation (23) and substituting the result 
into equation (22) yields 


D0, 


g 


3 g 


,, w Quantity plotted in fig. 7( a )> ref. 4 

*“ • X — ' — T 


1 + 3k 


.2 


(24) 


(l + k ' 2 ) 2 


This relationship is plotted in figure 2 for the various values of 0' 
used in reference 4 and as a function of reduced frequency u>' , where 


a)' = ^ = 0'k' 
U 


(25) 


In a like manner the relationship between the power spectra of 
yawing gusts and side gusts may be determined. From equation (ll) , 





~2* C l( u g) 

W 


(2 6) 


Dividing through by the power spectrum of the side gusts or vertical gusts 
as measured at the center of gravity gives 


*»g W 2 * C l( U g) 


K) 




w 


where at any given frequency 


I>|f c 


4U 2 ^ c t( u g) 

( C 0 W 2 


(27) 


Again, from reference 4, the power spectrum of the rolling -moment coeffi- 
cient due to horizontal gusts on the wing is given by 
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u * 2 l (%)w V 

^Cifu 1 = * Quantity plotted in fig. 9(a), ref. 4 

LK 6 JtU^ 

® Cz( k ') 

The quantity plotted in figure 9(a) of reference 4 is , 

u g^ L< "Zp a o 

where a 0 is trim angle of attack. Dividing the above equation by 

equation (23) and substituting the result into equation (27) yields the 
required relationship between the yawing- and side -gust spectra: 



4 x 


Quantity plotted in fig. 9(a), ref. 


1 + 3k' 2 
(l + k'2) 2 


4 


(28) 


This relationship is plotted in figure 3 as a function of cjd' for a 
range of values of 0 ' . 


Presentation of the data of figures 2 and 3 in the form of ratios 
between the gust spectra is not meant to imply that these data are 
independent of the approximate spectra chosen for 0g. These ratios 

represent the filtering effect of the wing on gusts having the frequency 
spectrum given by equation (23), and, if other approximations to this 
frequency spectrum were used, the filtering effect would not necessarily 
be the same. 


DISCUSSION 


The foregoing development leads to a set of equations which hold 
for small disturbances about some trimmed flight condition. In the appli- 
cation of the method, it is necessary to know the flight conditions, the 
stability derivatives of the airplane at those flight conditions, and some 
basic physical dimensions of the airframe. Such quantities are required 
for any dynamics study of an airplane and, aside from the equations and 
figures given herein, no additional information is required to obtain the 
lateral response of the airplane in power-spectral form to continuous 
atmospheric turbulence . 

As an aid in setting up the required calculations, a list of columns 
and steps which may be followed in the calculation of $0 is given in 

table I. The column numbers and headings axe listed vertically in this 
table; however, they would be arranged horizontally across the top of an 
actual calculation sheet. Although the calculation of is used as an 
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example, the spme steps are followed for the response of the airplane in 
any angular displacement, rate, or acceleration. For a given airplane and 
flight condition the first 12 steps (part (a)) may be tabulated from the 
equations and plots of this report. Steps 13 to 45 (part (b)) are used 
for the calculation of the response of the airplane in any lateral angular 
displacement (for the example, 0 is used). It is necessary to tabulate 
the real and imaginary parts of three transfer functions relating C j, 

C n , and Cy to that response parameter ^columns (l^) to (l|) J and to fol- 
low the indicated steps . For each of the other response motions or their 
derivatives it is only necessary to substitute the appropriate transfer 

functions into columns (l^) to (S) and to repeat the process . 

In calculating the frequency-dependent stability derivatives 
^columns (k) to care should be taken to insure that, when o> = 0, 

these derivatives agree with the steady-state aerodynamic derivatives 
used in calculating the transfer functions of equation ( 15 ). 

The method requires a choice of v alue s for the scale of turbulence L 

p 

and the mean square of the turbulence w g . Very little information is 
available at present on the proper magnitude for L; however, it appears 
to be within the range of 1,000 to 2,000 feet and probably closer to the 
lower figure . The mean square of the turbulence depends on the severi ty 

P 

of the turbulence to be considered. As an approximation, a value of w g 
of (3 ft/sec) 2 for light turbulence, (6 ft/sec) 2 for moderate turbulence, 
and (10 ft/sec) 2 for thunderstorms may be used. (See ref. 7*) 

In the present method the effects of turbulence of the airplane have 
been separated into the equivalent effects of rolling, yawing, and side- 
slip of the airplane in still air. The application of this concept, 
however, has been made in such a way that the effects of the u-, v-, and 
w-components of the gusts are taken into account as was done in refer- 
ence 1. Although the procedure of this report differs from that of 
reference 1, exactly the same effects are treated in both methods and 
both methods yield the same results. This agreement has been verified by 
using the present method to calculate responses for the example airplanes 
of reference 1. Since the responses of three airplanes are presented in 
reference 1, no numerical examples are presented herein. Those interested 
in the trends and relative effects of the different gust components over 
certain frequency ranges will find these effects discussed therein. 

One distinction should be noted between the present method and the 
more conventional treatment as contained in equation (3)- In the con- 
ventional treatment, yawing gusts are assumed to include both the rota- 
tional effects introduced by gradients of side gust along the fuselage 
and by gradients of horizontal gusts across the wing span. Thus, the 
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values of C 7 , C n , and C Y of equation (5) would be those for the 
t r Up 1 r 

entire airplane. In the present method, the yawing gusts are assumed 
to include only the effects of horizontal gusts on the wing, and the 
values of (^l r ) and (^n r ) used in equation (13) are those for the 

wing alone. Rotational effects of the side gusts along the fuselage 
are completely accounted for by the complex stability derivatives 
( c l ^ > ( C n \ > 811(1 ( C Y ) * 11118 method ls used because i-t allows 

V p/WT V P/FT V P/ft 

the utilization of the more accurate calculation of the fuselage pene- 
tration effects given in reference 8 and because the separation of the 
effects of the wing and fuselage allows the random distribution of hori- 
zontal gusts across the span to be taken into account. 

With most current airplane configurations in flight at low angles 
of attack a simplification of equations may be obtained by neglecting 
the yawing gusts on the wing. In the numerical examples of reference 1 
it was observed that this gust component had a negligible effect in all 
degrees of freedom for the airplanes investigated. This result is the 
basis for the assumption made in the derivation of the method used herein 
that the cross power between the gust-velocity components may be neglected. 
By the theory of isotropic turbulence only the cross power between the 
components herein referred to as yawing gusts and side gusts exists . How- 
ever, when the yawing -gust contribution to the motion of the airplane is 
small, neglecting this cross power appears to be justified. 

It was also found in reference 1 that the further simplification of 
neglecting the contribution of both rolling and yawing gusts in calcula- 
ting the response in sideslip is justified. 


The plots of the ratios of rolling gusts and yawing gusts to side 
gusts derived in this report give a physical picture of the relative 
iTrmnrtance of these eust disturbances at various frequencies. For cal 

*'1' a — 'Y rr 

culation purposes, the plots of the quantities and -r^ as 



2 

and 

rxr g 

0g 




functions of frequency on log-log paper (figs. 2 and 3) are convenient, 
but a plot of this type gives a somewhat distorted picture of the true 

variations of these quantities. For this reason, plots of 


_£ 


6 


and 



as functions of frequency on linear scales are given in figure 


Frequency is plotted in terms of the ratio of wing span to gust wave- 
length b/A. These curves show that, for the small values of 3' 
(ratios of wing span to scale of turbulence) ordinarily encountered. 


4 . 
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the ratio 




or 


0 g 


falls on a single curve as a function of 


b/A, except at low values of b/A. The curves reach a peak in the 
neighborhood of b/A = 0.5 to 1.0, as might be expected on the basis 
that the spanwise averaging effects would become important when the 
gust wavelength is shorter than the span. The hypothesis, expressed 
in reference 9 , that the gust gradient measured along the flight path 
would give an approximation to the effective spanwise gust gradient is 
indicated by the dashed line drawn on the figure (fig. 4(a)). The 
agreement in the trends of the curves with this dashed line shows 
that this hypothesis may be a reasonable explanation of the mechanism 
of the rolling and yawing effects of turbulence over an intermediate 
range of wavelengths. At shorter wavelengths, the spanwise averaging 
causes a decrease in the effects; whereas at long wavelengths, or low 
frequencies, another mechanism apparently comes into play to increase 
the rolling and yawing effects . This mechanism is believed to be the 
chance encounter of the wing with rolling and yawing gusts distributed 
along the flight path at relatively long intervals . This effect 
increases with the ratio of wing span to scale of turbulence. Possibly 
this result occurs because a wing of larger span effectively samples a 
larger portion of the atmosphere and is therefore more likely to 
encounter rolling and yawing gusts. This effect might be important in 
explaining increased lateral -control difficulties of large airplanes 
during landing approaches in rough air, if it can be shown that the 
scale of turbulence decreases at low altitudes. 


CONCLUDING REMARKS 


A procedure is presented for calculating in power-spectral form the 
lateral response of airplanes to random atmospheric turbulence. By fol- 
lowing the tabulated sample calculation procedure, these calculations 
may be made in a routine manner without detailed knowledge of the 
derivation of the method. It has been verified by using the present 
method to calculate the responses of the example airplanes of NACA 
Technical Note 3954 that the present method gives results identical 
to those of the method therein. The present method requires simpler 
calculations and provides a clearer physical picture of the relations 
between the various sources of lateral gust disturbances . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 15, 1957* 
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APPENDIX A 

FREQUENCY -DEPENDENT STABILITY DERIVATIVES 
DUE TO SIDE GUSTS 


All derivatives are referred to the center of gravity of the air- 
plane. The coefficient of rolling moment due to side gusts on the wing 
and vertical tail is given by 

A t co\ 

KHwt ■ N)w + (%)t K 1 + D* (A1) 

The coefficients of side force and yawing moment due to side gusts along 
the fuselage and vertical tail (see ref. 8) are expressed by 


p ( \ 

= 2n . 

2s 0 2 

/ x ik ol 

1 - 11 - ik A )e 

I^Yp^ 


c 
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0 

ro 

J. 1 -L ^ 


' s i - s o \ 2 
,k 2 - k l/ 
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’ -ikp -ik 2 

e - (l - ikp + ik 2 )e 


(A2) 




,e ik ° - 2i 


[ ko 3 

(xg - xp) (sp - Bp) 2 

(k2 - k i) 3 

-ikp 

(kp - 2i)e 


2k 2 - kp - 21 - ikpk 2 + I e 


2U‘ lk 2 


where 


k _ ^n 
kn ‘ "U~ 


(A3) 


(n=0,l,2) 


and where x Q , Xp, x 2 , s Q , and Sp are the profile dimensions of 
the fuselage and vertical tail. These profile dimensions, as given in 
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the derivation in reference 8, are illustrated in the following sketch: 


The fuselage is an ellipse truncated at the center of gravity, 2s q 
being the length of the minor axis and 2 xq being the length of the 
major axis. The vertical tail is a right triangle of base - x-^ 

and height - Sq. 

At a) = 0, equations (Al), (A2), and (A3) become, respectively, 


These quantities should be made numerically equal to the steady-state 
stability derivatives as obtained from flight-test measurements or more 
exact theories by suitably ad jus ting the dimensions of the assumed 
fuselage-tail profile. 




(A4) 



' ib[‘ 2so2x ° + (si ■ 8 ° )2 ( X2 + ?)] (A6> 
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APPENDIX B 

AIRFRAME TRANSFER FUNCTIONS 


In coefficient form, the lateral moments and forces on an airplane 
are related to the angles which define the lateral motion by 




0 


c i 

< \|r 

► = < 


3 




where £a[] is defined in equation (4) for the case where the X-axis 

is initially alined with the relative wind (stability axes). The change 
in the angles due to changes in the moment and force coefficients is 
then expressed as 




Cz 

< \|f 

II 

1 1 

1 

H 

0 

3 

0 


Cy 

V. 


One method of obtaining the inverse of a matrix |^aJ is given by the 
relationship 



Adjoint [A T j 


where the prime denotes the transpose of [a] . It may be seen that, 
when [a] is defined by equation (4), ^a] - - 1 - will be made up of 

transfer functions as given in equation (15). The transfer functions 
are defined by 

1 + C nr )D 2 + (|c nr C Y|3 + 2\iC n ^ - |c n(J C Yr )D - 


c 


n 0 


tan 7 
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V 2 K xz d5 + - 2 K X Z c Y 0 ) I)2 + 2 ( C n 0 C Y p - c Y 3 C np) D + C L c n 3 

= |M- _2K XZ ^ “ ^Yj.) + ^^Ypj 05 ' " 2 ^ C l( K XZ tan 7 + K Z 2 ) + 
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A | = 8(1 3 ( k x \ z 2 - - 2(i 2 |K x 2 ^2K z 2 C Y p + + Kg 2 ^ + 
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APPENDIX C 

TWO EQUIVALENT WAYS OF EXPRESSING THE POWER SPECTRAL 
RELATIONSHIP BETWEEN GUSTS AND 
AIRPLANE MOTIONS 


The statement is made in the main text that, in the derivation of 
the power spectral relationships between the airplane motions and the 
gust-velocity inputs, the same result is obtained whether the matrix 
relationships are treated in the form 

A _1 (a>)| 2 |g(o))| 2 (Cl) 


or in the form 



(C2) 


In the former case, cross -power terms between the elements of the 
matrices (i.e., between the moments and forces on the airplane due to 
each gust component acting on the various parts of the airplane simul- 
taneously) will appear, whereas in the latter case, these quantities 
do not appear explicitly but are taken into account when the matrices 
are first multiplied together. The former approach is used in the 
method of reference lj the latter and simpler approach is used in the 
method of this paper. In order to show the equivalence of the two 
approaches it is necessary only to prove equality between equations (Cl) 
and (C2), that is, that the square of the absolute value of the product 
of two matrices is equal to the product of the squares of their absolute 
values . Since matrix equations may be treated as linear algebraic quan- 
tities, the matrix equality defined in the main text 


[ a(oo)] _1 [5(a>)] = [a _1 G(cd)] 

will be denoted by the complex quantities 

AB = C 

The complex conjugate of this equality may be shown to be 


(C3) 

(C4) 

(C5) 


A*B* = C* 
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Multiplying equations (C4) and (C5) yields 

AA*EB* = CC* 


which is equivalent to 


A| 2 |Bp . |c| £ 


Hence , 


-1~, x|2 


| A -1 (o>) | | G(o>) | = |a -1 G(cd) 


(C6) 


where the elements of these matrices are likewise complex quantities . 

Insight into the difference between the two approaches may be shown 

by considering one element of the [a _ 1 gJ matrix. Expanding equation (C3) 
by means of equations (13) and (15) of the main text gives 


ill! 


2 (^p) w 


'n 


;J j !(%)«’ 


, A. 


or, in its expanded form, 


D0 g 2\ P/W Cj 2 \ n P/ W C n 


A+ll 


A 


(CT) 


If the real and imaginary parts of equation (C7) are grouped, then the 
absolute value squared of equation (C7) becomes 


D 0g 


R 


1 (c 7 ^ A + i/c \ -^-ll + -l[ifc l \ A + i/c ' ) A 

2 V l P/W C z 2 \ "pjw C n J [_2\ l pjw C z 2 \ Hp/ W 


(C8) 


A different form involving cross -power terms may be obtained by 
writing the conjugate of equation (C7) 
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0 1 




0 


D0j 


g 


and obtaining the product of equations (C 7 ) and (C9): 

4-(4T = rf c i 1 2 + HO 2 

D0g\D0g / H*P/W C^cJ 4V“pJw C n lc n 


(C9) 


Since 
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equation (CIO) may be written in the form 
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D0, 
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■ l( c 'p)w £ 
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+ 2R 
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The term 2R 


*(%)«(%)« £(£) 


(cn) 

is the cross -power term between the 


rolling and yawing moments on the wing, and had the coefficient (c^- \ 

\ *P/W 

not been zero two other cross -power terms would have appeared. 

The expressions of equations (C8) and (Cll) are equivalent. Either 
form may be used, but, for the purposes of this report, equation (C8) 
is more useful and the illustrative tabulation process has been set up 
on this basis. 
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TABLE I.- SAMPLE CALCULATION PROCEDURE 
(a) Calculation of frequency -dependent derivatives and gust spectra 


Column 

number 


Column heading 


Instructions 
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© 

© 

© 

© 

© 

© 

© 
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03 
CD 1 


(%)wt] 
( Cl p) wtJ 
[( Cn p) ft] 

( Cn 3) FT_| 
|( Cy 3) FT | 
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G) — 

u 


= <d'/0' 


From equation (Al) 


From equation (A3) 


From equation (A2) 


Quantity plotted in fig. 2 


Quantity plotted in fig. 3 x 


Quantity plotted in fig. 1 x 


'a£l v \ 2 
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(b) Calculation of roll response 


Column 

number 

Column heading 

Instructions 

© 

R(0/Ci) 
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© 
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Evaluation of airframe transfer functions 
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* (appendix B) at specified values of cd' 
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TABLE I.- SAMPLE CALCULATION PROCEDURE - Concluded 
(b) Calculation of roll response - Concluded 
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number 

Column heading 

Instructions 
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Figure 1.- Power spectral density of the side-gust component of isotropic turbulence 
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Figure 2.- Ratio of rolling-gust power spectrum to side-gust power spectrum as a function of 
reduced frequency for various values of (3 ' , the ratio of wing span to scale of 

turbulence . 
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Figure 3 .- Curves for determining the ratio of yawing -gust power spectrum to side-gust power 
spectrum for various values of the ratio of wing span to scale of turbulence (3 , wing angle 
of attack c^, and wing stability derivatives C and C Zp . 
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A 

(a) Ratio of the rolling -gust component to the side-gust component. The approximation based on 
a constant antisymmetric gust-velocity gradient over the wing span is indicated by the 
dashed line. 


Figure k.- Ratio of the absolute amplitudes of the rolling- and yawing -gust components to the 
side-gust component as a function of the ratio of wing span to gust wavelength. 
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(b) Ratio of the yawing -gust component to the side-gust component. 

Figure k.- Concluded. 
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